On frequency errors of nanomechanical-resonators-based-on quantum computing 
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We study the consequence of the frequency errors of individual oscillators on the scalability 
of quantum computing based on nanomechanical resonators. We show the fidelity change of the 
quantum operation due to the frequency shifts numerically. We present a method to perfectly 
compensate for these negative effects. Our method is robust to whatever large frequency errors. 
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I. INTRODUCTION 

Scalability is one of the most important issues in the re- 
alization of quantum computing. For example, in order 
to factorize a 200-digit number using a quantum com- 
puter, one needs to manipulate thousands of qubits 
However, in realistic physical systems, there are many 
imperfections, such as the quantum decoherence, the de- 
vice errors, and so on. These imperfections seriously limit 
the power of a quantum computing device. In particular, 
in a large scale quantum computing, the small error of 
each individual device can accumulate and this may lead 
to the failure of the final result Q . 

In the recent years, methods for scalable quantum com- 
puting based on artificial quantum systems have been ex- 
tensively studied A promising scalable quantum 
computing architecture based on spin system of nonomc- 
chanical resonators (NAMRs) was proposed by P. Rabl 
et. al. 0. The spins have a long decoherence time 
and the NAMRs can be fabricated on a large scale. The 
quantum motions of the NAMRs can strongly interact 
with the spins J7[ and induce strong couplings between 
the spins @, S [3 • 

However, the inevitable fact is that the frequencies of 
the NAMRs can't be exactly the same in practice. There 
are always frequency errors due to imperfect fabrications. 
Experimental results show that there can be ±1.0% de- 
viations from the averaged frequency [loj . These quan- 
tum imperfections may cause exponential suppression of 
quantum computations |ll| , thus it is an important issue 
to find out quantitatively the impacts of these frequency 
errors and the efficient method to compensate for the 
negative effects due to the device errors. 

In this article, we consider a scalable quantum comput- 
ing architecture consisting of TV spin qubits whose inter- 
actions are mediated by an array of N NAMRs [1, . The 
frequency errors of the NAMRs shift the frequencies of 
the collective modes, the NAMR-spin coupling strengths 



'Electronic address: |xbwang@mail. tsinghua. edu.cn ] 



and result in fluctuations in the spin-spin couplings. For 
a given evolution time, the spins' final state depends on 
the coupling strengths among them. Supposing a cer- 
tain duration is needed as the right evolution time to 
produce a certain state in the ideal case of no frequency 
errors, the same evolution will produce a wrong state in 
the actual case with frequency errors, due to the fluctu- 
ations of the spin-spin couplings which would reduce the 
quality of the quantum computing in the spin-spin inter- 
actions mediated by the NAMRs. Here, we analyze the 
impacts of the frequency errors on the quantum operation 
fidelity and figure out a limitation of the scalability of the 
quantum computing architecture. Our numerical simula- 
tions show that the quantum operation fidelity decreases 
rapidly with the number of the NAMR-spin elements in 
the present of frequency errors and these frequency errors 
finally limit the scale of the system. Based on the analy- 
sis, we propose a method to compensate for the negative 
effects of whatever large frequency differences. 

This article is arranged as follows: In Sec. |ll]the model 
is given first and then the method of diagonalizing the 
quadratic boson Hamiltonian is briefly reviewed. Next, 
in Sec. IIIIl we show the consequence of frequency er- 
rors on quantum operation fidelity by numerical simula- 
tion. We study how the influence of the frequency errors 
change with the number of the NAMRs TV. In Sec. IIVI 
we present our theoretical method to resist the frequency 
errors. Finally, we present discussions and conclusions in 
Sec. El 



II. MODEL AND METHOD 
A. Model 

As shown in Fig. [1] we consider a system consisting of 
an array of iV NAMRs, which are charged and interact 
capacitively with nearby wires interconnecting them @ . 
A magnetic tip is attached on the free end of each NAMR. 
An NAMR with fundamental frequency Wi and effective 
mass m magnetically couples to an electronic spin qubit 
associated with a nitrogen- vacancy (N-V) center located 
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FIG. 1: (Color online.) An NAMRs-based spin chain. The 
NAMRs vibrate in the direction. Each two nearest- 
neighbour NAMRs are connected by an individually iso- 
lated wire and interacted capacitively with each other. Each 
NAMR is attached with a magnetic tip on the free end and 
magnetically couples to a spin qubit (N-V center) located in 
the substrate below, which is driven by a local microwave. 



between the spin i and the spin j, which could be used 
for scalable quantum computation @ . 



B. Diagonalizing the quadratic boson Hamiltonian 

When the frequencies of the NAMRs are exactly 
the same (denoted by Wr) and the nearest-neighbour 
NAMR coupling is a constant g, the frequencies of col- 
lective modes can be obtained by solving the eigen- 
value equations analytically and are given by u;„ = 
y^uj^ + 4u}ig{l + cos[(n -f l)7r]/A^}. In practice, the fre- 
quencies of all the NAMRs can't be exactly the same due 
to imperfect fabrications. These systematic frequency 
errors of each NAMR are independent and are not the 
same, the frequencies of the collective modes of N cou- 
pled NAMRs u!n can't be simply calculated by the pertur- 
bation method. A useful method is described as follows. 

The Hamiltonian ([T]) can be rewritten as 



in the substrate below [7|, |12| . Each spin is driven by a 
local microwave to form a pair of dressed states in or- 
der to match the NAMR frequency 01 . The interaction 
Hamiltonian between the NAMR and the corresponding 
spin qubit is iJ*,. = ^{a] + ai)a\ Here a\ (ui) is 
the natural creation (annihilation) operator of the fun- 
damental vibrational mode of the i*'' NAMR. al is the 
Pauli-z operator of the i*^ spin. The coupling strength 
A = gsfJ-BGin0.a/fi with gg = 2, the Bohr magneton /iB, 
the magnetic field gradient Gm and the amplitude of zero- 
point fluctuations oo = ^y h/{2muji). 

The Hamiltonian of these N coupled NAMRs is (set- 
ting /i = 1) [I: 

N ^ N 

with gij = al[d'^Uci/{dzidzj)]\zi=o,z,=o- Here C/d = 
-(J72/2)CsCw/(Ce + Cw) is the electrostatic energy be- 
tween two nearest-neighbour NAMRs connected by a 
wire of self-capacitance Cw with Cs = Ci{zi) + Cj{zj), 
Ci{zi) ~ C(l — Zi/h). C/v is the applied gate voltage on 
each NAMR, h is the electrode spacing, Zi is the tip po- 
sition of the NAMR and C is a constant. The first item 
of iJph is the free Hamiltonian of the NAMRs and the 
second item is the interactions between them. The in- 
teraction gi j consists of two parts: the self-coupling gi^i 
and the coupling between different NAMRs gi j [i ^ j). 

If the Hamiltonian iJph is diagonalized by defining col- 
lective modes, the total NAMR-spin coupling Hamilto- 
nian (the sum of Hl^. over i) can be rewritten as the cou- 
pling between the z-components of the electronic spins 
and the collective modes of the NAMRs. After a trans- 
formation, the effective spin-spin interactions mediated 
by the NAMRs can be obtained and have the form 
^^off = Tlii j ■^i-j'^\'^z with Mi,j the coupling strength 



N 

Hph = a^S^a = Y,^kbibk, (2) 
fe=i 

with 6^, (6fc) the creation (annihilation) operators 
for the phonons of the collective modes and a = 
[oi, • • • , flAT, a|, • • • , ajy]-^ the natural creation (annihila- 

\ A B 

tion) operators of all the NAMRs. Here ^ = ^ ^ ' 

A is an X iV diagonal real matrix depending on the 
free Hamiltonian of these NAMRs and i? is an A'' x 
symmetric real matrix governed by the coupling model 
between NAMRs. 

Introducing an auxiliary matrix D = — B^, the 
para-values and para-vectors iTsI (and also the diag- 
onalizing para- unitary matrix |14l |) of ^ can be con- 
structed from the eigenvalues (denoted by with 
k ~ 1---N) and eigenvectors (denoted by ^k) of 
D [Hi. Then the 2A^ x 2A^ transformation ma- 
trix is given by = [C1C2 ■ ■ • CnQn+i ■ ■ ■ Gat] with 
Cfc and C,N+k [lal are 2A^-column vectors The 
transformation matrix should satisfy {,P)~'^^^J~^ = 

J- jl/2 .1/2 ,1/2 ,1/2 ,l/2\ , 

ai&g[di , • • • , , d]^ 1 ■ ■ ■ i ) preserve 

the commutation relations of bosonic operators. The col- 
lective frequencies of these NAMRs are 

i:^k = 2d]!\ (3) 

By denoting b = • • • , bj<4,b\, • • • , fojy]^, we obtain the 
relation a = J~^h. The collective frequencies are shifted 
by the frequency errors and the numerically result for 
A^ = 11 is as shown in Fig. [51 

For the system we considered, we only need to consider 
the nearest-neighbour NAMR-NAMR interactions. The 
coupling strengths between nearest-neighbour NAMRs 
are assumed to be the same and denoted by gt^i+i = g. 
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FIG. 2: (Color online.) The collective frequencies of the 
NAMRs are shifted by the random frequency errors for — 
11 with numerical calculation. The parameters used here are: 
Wr/(27r) = 1 MHz, the coupling strength: fl/(27r) = 500 kHz. 
The left panel (a) is the result by running the program one 
time while the right panel (b) is the average by running the 
program 200 times. The offsets of the collective frequencies 
are optimistic estimations. 
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The matrix A and B are given by 
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FIG. 3: (Color online.) Fidelity decreases as the number of 
the NAMRs increases by numerical simulation (we run the 
simulation program 200 times for average to get the tendency 
of the fidelity with the NAMR number). The memory and 
time needed for simulating quantum systems exponentially 
increase as the scale of the quantum system N increases. We 
only simulate this system from N = 2 to N = 11 on an 
ordinary personal computer. The frequency errors compared 
with the central frequency are: 1.0% (black solid square), i.e.. 
Am = 1.0% and the same below, 3.0% (blue solid circle), and 
5.0% (red solid star). 



III. FIDELITY ANALYSIS 

Based on the model and method above, we analysize 
the effects of the frequency errors. We shaU calculate the 
fidelity of the assumed state from the ideal model without 
any frequency error and the actual state from the more 
realistic model with independent frequency errors in each 
individual NAMR. 

Denoting the frequency of each NAMR = uj^ + ei 
with Wr the averaged frequencies of all the NAMRs and 
gj the frequency error, the coupling strength between 
the NAMR and the corresponding spin can be written 
as Aj = Ai/wr/wi, which means that the frequency er- 
rors generally fluctuate the NAMR-spin couplings. The 
NAMR-spin coupling can be rewritten as 

N . N 

Y^^(a., + al)al= ^ A„..(6„ + fot )a:, (5) 

i n,i—l 

with 

An,i — ^i{Ji^n + Ji,n+N)^ (6) 

the coupling strength between the n"^ collective mode 
and the z component of the i"^ spin. Here {>/i~/} are 



the matrix elements of J ^ . Therefore the effective spin- 
spin coupling is M^j = Y,n'^n,i>^n,] / {^^^n)- For a given 
evolution time tg, the spin-entangling operation can be 
described by Ug{tg) = exp{i{J2i Mijalai)tg} 

For N NAMRs with an initial state \tpin) = [(|0) + 
|l))/\/2]®^, the relative frequency errors are denoted by 
Ai = £i/wr. After a given evolution time tg, the spins' 
sate becomes \(j)o{N,{Ai})) = Ug{tg)\ipin) ■ The fidelity 
between the actual state and the ideal state is 

F{N, {A,}) = \{MN, {A,})\MN, {A, = 0}))|. (7) 

We consider an array of Si-NAMRs with an averaged 
frequency a;r/(27r) = 1 MHz, oq « 1.86 x 10^^'^ m and 
g/{2n) 500 kHz ^. For gradient Gm = 9.6x 10*^ Tm"!, 
the resulting coupling strength between an NAMR with 
frequency uj^ and a spin is about A/(27r) w 50 kHz. In this 
case, the effective nearest-neighbour spin-spin coupling 
is about a few kilohertz. We choose the evolution time 
tg = 0.3 ms to complete one quantum operation. (For 
example, given = 2, if we want obtain the entangled 
state (|00) + |11) +i|01) +z|10))/2 from the initial state 
[(|0) + |1))/a/2]®^, the evolution time and the two spin 
couphng strength M should satisfy = 7r/4|M| [1]. The 
spin dccoherence time r2 ~ 6 ms [13 and the dephasing 
time induced by nuclear-spin fluctuations T2 ~ 0.35 ms 
for N-V centers [l^ are observed experimentally. ) 
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FIG. 4: (Color online.) Fidelity F is approximately linear 
with the number of the NAMRs A*". The frequency errors 
here satisfy Ei/ur G [—1.0%, 1.0%]. The black dots are the 
simulation data and the red line is the fitted curve. 



Our numerical simulation shows that the fidelity 
F{N,{/^i)) decreases with the number of the NAMRs 
N, as sketched in Fig. [31 In Fig. |31 the red solid 
star gives F(2,Am = 5.0%) w 0.988 while F(11,A„, = 
5.0%) « 0.971. F(2,A,„ = 1.0%) w 0.9996 while 
i^(2,Am = 5.0%) w 0.9879. Here F{N,A,^ = 5.0%) 
is the value for F{N, {A^}) where the values of each A^ 
are randomly chosen from the range [— Ami+Am]. This 
shows that, as the number of qubits increases, the errors 
in the target state rise, hence one may end up with a 
wrong result with large probability in a large scale quan- 
tum computation, even though the frequency errors of 
each individual NAMR are small. 

Fig. [3] also shows that F is approximately linear with 
A''. For example, suppose A„i = 1.0%. For A^ = 2 • • • 11, 
the linear fitting gives the relationship between F and N: 
F = -0.000080952 x N + 0.99976 with the correlation 
coefficient r = 0.99 as shown in Fig. For a given Fq, 
if wc want the fidelity F > Fq, the length of a spin chain 
should be less than about [10^ x (0.99976 -Fo)/0.80952j 
spins. The symbol [P\ denotes the maximal integer not 
great than P. For example, wc obtain iV < 4.6 x 10^ for 
i^o = 2/3. It can be seen that the frequency errors do 
limit the maximum length of the spin chain if we want 
to produce a faithful resultant state and finally limit the 
extensibility of the fault-tolerant quantum computing. 



IV. COMPENSATION METHOD 

As analyzed above, the frequency errors can lead to 
fidelity decease of the quantum state evolution. The de- 
viation cannot be compensated by simply adjusting the 
evolution time tg when A^ > 2, as shown in Fig. [5l In 
what follows, we propose a method to solve this prob- 
lem. The main idea is that the errors caused by fre- 
quency differences can be perfectly compensated for by 
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FIG. 5: (Color online.) The fidelity F changes as a function 
of the evolution time in the case of frequency errors for a 
fixed A. We set the system state at the time = tgo = 0.3 
ms without frequency errors as the ideal final state. F denotes 
the fidelity between the ideal final state and the actual state 
at an arbitrary time tg. The infidelity (1 — _F) caused by fre- 
quency errors can't be compensated for by only adjusting the 
evolution time when A > 2. In this figure, the frequency 
errors are within 1.0% compared with lo^. The black (solid) 
curve is for N = 2, the dashed (red) curve is for N — 5 and the 
dotted (blue) curve is for A = 8. The curves in panel (b) are 
the plot details of panel (a) for the regime tg £ [0.290, 0.315] 
ms. 



controlling the interaction time of each two adjacent spin 
qubits block. With this compensation method, one can 
produce arbitrarily large-scale states with whatever large 
frequency error of each NAMR. 



For a one-dimensional N NAMR-spin chain, an auxil- 
iary switchable control voltage gate is added between ev- 
ery two nearest-neighbour NAMRs to obtain switchable 
coupling between them, as depicted in Fig. |6l When the 
control voltage is zero, the nearest-neighbour NAMRs 
couple with each other and induce effective spin-spin in- 
teractions; and when the control voltage is switched on 
and equals to the NAMR's gate voltage, the coupling of 
the two NAMRs is switched off, hence the NAMRs induce 
no effective spin-spin interactions. In addition, we as- 
sume that we have obtained the fundamental frequencies 
of the NAMRs by average after multiple measurements 
(the systematic errors still exist due to the limitations of 
the fabrication techniques and the random measurement 
frequency errors arc eliminated by average) . The effective 
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FIG. 6: (Color online.) Switchable coupling between two 
NAMRs (top view). A switchable control gate (dashed box) 
is added between every two nearest-neighbour NAMRs. For 
a control voltage Uc ~ 0, the two NAMRs couple with each 
other. For Uc = Uv, there is no energy associated with charge 
flowing from the gate capacitors onto the wire and thereby 
the coupling between the two NAMRs is switched olT. 




FIG. 7; (Color online.) A sketch for the procedures to com- 
pensate for the frequency errors. 



couphng strength between every two nearest-neighbour 
spins in the case of frequency errors can be theoretically 
calculated by the formula M^^i+i by assuming that the 
other NAMRs and spins do not exist. 

For simplicity, we assume N is an even number (for 
an odd N, we only need to regard the last spin as an al- 
ready entangled spin pair). One can achieve the perfect 
result through the following procedures as sketched in 
Fig. [T] (1) switch on the coupling between the (2i — 1)*'^ 
and the (2i)*'^ spins for z = 1 • • • N/2 and switch off the 
coupling between the (2i)*^ and the (2i 4- 1)'^ spins for 
i = 1- ■ ■ {N/2 — 1). (2) control the evolution times be- 
tween each two-spin pair to obtain the ideal entanglement 
between them and then switch off the coupling, respec- 
tively. (3) switch on the coupling between the {2iY^ and 
the (2i -I- 1)*'^ spins for i = 1 • • • (A^/2 - 1) when all the 
couplings are switched off in (2). (4) repeat (2). The to- 
tal time needed for all these processes is approximately 
two times of that in the case of no frequency errors. In 
addition, this method also avoid the influence of effective 
sub nearest-neighbour interactions between the spins. 

Proof. The operators of different spins communicate with 
each other. In the case of no frequency errors, the 
NAMR-spin couplings reduce to a constant and denoted 
by A. All the effective nearest-neighbour spin-spin cou- 
plings induced by the NAMRs are equal to each other 
and denoted by M. According to the Ising model, the 
evolution operator of the spin chain for a given evolution 
time tg is: 



N-l 



(8) 



= Yl exp{iMcr>*+4g} 



In the case of frequency errors, the NAMR-spin cou- 
plings and the collective modes of the NAMRs are shifted 



by the frequency errors. Denoting the coupling strength 
between the z*^ and {i + 1)**^ spins as Mi^i+i, the ideal 
entanglement between two spins can be obtained by ad- 
justing the evolution time in the two spin case as shown 
in Fig. [5]and the time needed is denoted by When 
all the four steps are completed, the evolution operator 
of these spins can be described by: 



N-l 



1,2 •• • 



1,N 



) = Yl exp {iM,,,+iaicTl+H,,,+i} 



(9) 



By precisely controlling the evolution time i^^i+i = 
(Mtg)/Mi^i+i, the evolution operators in the two cases 
are equal to each other and the two final states must 
be the same if the system evolutes from the same initial 
state. □ 



We should point out that the entanglement of the spins 
prepared by the following steps is not perfect: (1) switch 
on all the nearest-neighbour NAMR-NAMR couplings. 
(2) switch off the coupling when the nearest-neighbour 
spins involute to the maximum entanglement sates (the 
effective interaction between the nearest-neighbour spins 
can be numerically calculated), respectively. A spin 
interacts with all the other spins through the collec- 
tive modes of the NAMRs when the control voltage is 
zero. These interactions inevitably include sub nearest- 
neighbour interactions and so on. These sub nearest- 
neighbour interactions decay as Mi^i±m ~ (g/wr)*-™"^-* 
0, but their influences are generally comparable with 
those due to frequency errors. Unwanted entanglement 
between the spins which are not nearest-neighbours is 
generated and the Ising interaction produced by the 
above two steps is unsatisfactory. 
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V. DISCUSSIONS AND CONCLUSIONS 

Although only the one-dimensional case is discussed in 
the context, the similar analysis on the two-dimensional 
NAMRs-spins quantum computing architecture Q can 
also be handled by using the same method. In the two- 
dimensional configuration, NAMRs are ordered on a two- 
dimensional lattice and couple to their four neighbours 
electrostatically. In the two-dimensional case, the effects 
of the frequency errors can be analysed by adjusting the 
matrix B in equation ()4bp to be a general matrix instead 
of a tridiagonal matrix when considering the nearest- 
neighbour couplings. In the compensation method, a 
switchable voltage gate should be added between each 
two nearest-neighbour NAMRs. 



In summary, we study the impacts of the frequency dif- 
ferences on the scalability of a promising NAMRs-based- 
on quantum computing architecture. The influences on 
the quantum operation fidelity are analyzed in detail, and 
a method is given to compensate for the negative effects 
of these frequency differences. 
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